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ABSTRACT 

We explore the ratio (C/M) of carbon-rich to oxygen-rich thermally pulsing asymptotic giant branch 
(TP-AGB) stars in the disk of M31 using a combination of moderate-resolution optical spectroscopy 
from the Spectroscopic Landscape of Andromeda’s Stellar Halo (SPLASH) survey and six-filter Hub¬ 
ble Space Telescope photometry from the Panchromatic Hubble Andromeda Treasury (PHAT) survey. 
Carbon stars were identified spectroscopically. Oxygen-rich M-stars were identified using three differ¬ 
ent photometric definitions designed to mimic, and thus evaluate, selection techniques common in the 
literature. We calculate the C/M ratio as a function of galactocentric radius, present-day gas-phase 
oxygen abundance, stellar metallicity, age (via proxy defined as the ratio of TP-AGB stars to red 
giant branch, RGB, stars), and mean star formation rate over the last 400 Myr. We find statistically 
significant correlations between log(C/M) and all parameters. These trends are consistent across dif¬ 
ferent M-star selection methods, though the fiducial values change. Of particular note is our observed 
relationship between log(C/M) and stellar metallicity, which is fully consistent with the trend seen 
across Local Group satellite galaxies. The fact that this trend persists in stellar populations with very 
different star formation histories indicates that the C/M ratio is governed by stellar properties alone. 


1. INTRODUCTION 


Asymptotic giant branch (AGB) stars are important 
for understanding galaxies’ integrated light and resolved 
stellar populations. They are major contributors of near- 
infrared (NIR) flux, contributing ^ 20% of a galaxy’s 
NIR lig ht in the local universe, and up to 70% at high- 


redshift ([Boyer et al.|2011[ Melbourne et al. 12012 


2013 IVIelbourne Boyer|'2013" Villaume et al.'2015 


Conro’ 


f 


They also remain one of the least understood phases of 
stellar evolution, with outstanding questions regarding 
calibration of the thermally-pulsating AGB (T P-AGB) 
phase, dred ge-up, opacities and mass loss (e.g. jMarigo 
et al. |2013 , and references therein). 

AGB stars are broadly characterized by whether their 
atmospheres contain excess carbon or excess oxygen. 
Stars with free oxygen (C/0 < 1) are deemed M-type 
AGB stars, or M-stars. These stars undergo the third 
dredge-up (3DU), which brings newly formed carbon and 
s-process elements to the surface. This pollution alters 
the atmospheric chemistry, and causes the star to transi¬ 
tion first to S-type (C/0^ 1) and finally to C-type (C/0 
> 1). The C-type stars with free carbon are also known 
as carbon stars. 
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Because the transition from M-star to C-star depends 
on metallicity and stellar mass, the ratio (C/M) of 
carbon-rich to oxygen-rich AGB stars is a useful tool 
for studying the evolution of TP-AGB stars, and the 
galactic environment in which the stars formed. For ex- 


Local Group (e.g. Cioni & 

Habing|2003l |Battinelli & De- 

mers|2004bl 

Cioni et al.|2UU8l|Battinelli & Demers|2uu9al 

Boyer et al. 

20131), have been used to constrain models 

of AGB stars (iKarakas 2014) and metallicity gradients 

of the host galaxies (e.g. 

Gioni et al. 

20081 |Feast et al. 


2010 ). 

With its complicated star formation history (SFH) and 
metal-rich environment, M31 is a particularly powerful 
laboratory for continuing this work. Due to its size and 
distance, past work on the AGB population in M31 has 


(Brewer et al. 1995, 

1996 

Battinelli & Demers 2005), 

Hubble Space Telescope (H 

ST’) pencil beams (|Stephens 


et al.|2003||B oyer et al.|^2013 
observation ( |Davidge 200l| 


2013), or adaptive-optics assisted 


Davidge et al. 2005). For¬ 
tunately, recent large-scale surveys have opened up the 
AGB stars of M31 to the same scrutiny as elsewhere in 
the Local Group. 

In this work, we use data from both t he Panchromatic 
Hubb le Andromeda Treasury (PHAT; |Dalcanton et al. 
2012) and the Spectroscopic and Photome tric Landscape 
of Androme da’s Stellar Halo (SPLASH; |Guhathakurta 
et al. 112006) surveys. The pairing of HST' photometry 


and Keck spectroscopy across the disk of M31 allows us 
to not only distinguish carbon-rich and oxygen-rich AGB 
stars, but also to evaluate the C/M ratio as a function of 
a wide variety of parameters. We can, for the first time, 
self-consistently evaluate the C/M ratio across a variety 
of age, metallicity and SFH environments. 

Section describes our data, including the surveys 
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from which we obtain our spectra and photometry, and 
our method of distinguishing between carbon-rich and 
oxygen-rich AGB stars. Section examines the effect of 
differen^ossible AGB classifications on the G/M ratio. 
Section 0 presents the resulting C/M ratios as a function 
of radius, metallicity, and star formation history (SFH). 
Section discusses these results and possible biases. 

2. DATA 


2.1. Photometry 


The photometry used in this work is from the PHAT 
survey, w hich resolved ^ 117 m illion stars in the disk 
of M31 d Williams et al.| 2014|). Images were taken 
with UV (T'27514/ and Ta 3614/), optical (F475IT and 
F814IT) and NIR (FllOW and F160IT) filters, using 
HST WFC3/UVIS, ACS/WFC and WFC3/IR, respec- 
tively. The pho tometry was performed using DOLPHOT 
(Dolphin][2002), and “good star” (gst) cuts were made 
using the signal-to-noise (S/N), sharpness, and crowd¬ 
ing parameters output by the photome try pipeline. For 


furthe r det ails, we refer the read er to [Dalcanton et al. 
(2012) and Williams et al. (2014). 

When c ompar ing the C/lVl ratio to t he st ellar metal¬ 
licity in § 4.2.2 and the age pro xy in ^ a. we use the 
single camera ACS photometry (|Dalcanton et al.||2012p 
further correc ted for bias, comp l etenes s, and foreground 
extinction by Gregersen et al. (2015|) (hereafter G15). 
Photometric bias is caused by the effect of unresolved 


and bright stars in crowded regions, which leads stars 
to appear brighter and with a color closer to the mean 
color of the region. Completeness, the fraction of stars 
observed vs. the number of stars present, has a color 
dependent effect on the stellar population. While com¬ 
pleteness does not impact the magnitude of a star, it does 
impact quantities derived from the population (e.g. me¬ 
dian metallicity). To correct for foreground extinction, 
G15 assume Ay = 0.17 mag. For further details on these 
corrections, we refer the reader to the original paper. 


2.2. Spectroscopy 


mented elsewhere (e.g. Guhathakurta et al.|2006| 

Gilbert 

et al. 1120091 IDorman et al.|2012| [Gilbert et al.|20J 

.21 |Dor- 

man et al.|2015), but those specihe to the disk an 
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summarized here. 

Our spectroscopic targets were selected from optical 
photometry. The majority 63%) were targeted using 
HST photometry from PHAT. The rest were targeted 


using CFHT i' photometry (discussed in Gilbert et al 


20121). In all cases, stars were sel ected to be isolated, 


without close, bright neighbors (see Dorman et al.||2012 
for details). 

The disk dataset contains 10,61 9 optical spectra ta ken 
with the DEIMOS spectrograph (|Faber et al.||2003|) on 
the Keck H 10-m telescope. Approximately half (5o23) 
of these spectra were observed with the 1200 line mm“^ 
grating, which has a dispersion of 0.33A pixeR^ and a 
central wavelength of 7760A. The rest (5296) were taken 
with the 600 line mm“^ grating, which has a dispersion 
of 0.65A pixel”^ and a central wavelength of 7000A. 

The spectra were reduced using the spec 2d and spec Id 
software modified for the SP LASH survey (Gooper et al. 
2012 Newman et al. 2013). The spec2d routine firsts 



Spectral Type 


Fig. 1. — Distribution of spectral types represented in our suite of 
high S/N Milky Way template spectra. “C” denotes carbon stars. 


extracts a one-dimensional spectrum from the two- 
dimensional spectral data, and the spec Id routine deter¬ 
mines the redshift of this spectrum by cross-correlating 
with a series of templates. This cross-correlation also 
serves to distinguish galaxies and failed observations that 
returned no flux. 

We refer to these 10619 spectra as our “science spec¬ 
tra,” and they have a median S/N of 3.4 pix”^. 


2.3. Spectroscopic Templates 


To complement the Milky Way spectroscopic templates 
included in the spec2d pipeline, we observed a series of 
Milky Way carbon stars to use as radial velocity and 
spectral type templates. The stars were sele cted from 
the SIMbAd database (Wenger et al. 2000) based on 
magnitude, position, and spectral type. T'heir apparent 
magnitudes are in the range 7.74 > V > 6.21 and their 
spectral types span G5 to G8 based on the Keenan an d 
Morgan classification scheme (Keenan & Morgan] 1941). 

The observations were carried out using the 600 
line mm”^ DEIMOS grating centered at 7000A and the 
GG455 filter to block shorter wavelength light. We 
use the long-slit mask LVMslits with a 0.8” wide slit. 
This instrumental configuration provides a waveiength 


coverage of 4800-9500A with a spectral resolution of 
0.65A pixel”^. The raw two-dimensional spectra were 
reduced and extracted in the same manner as the sci¬ 
ence spectra. 

Our final suite of spectroscopic templates contains 31 
stars with a median S/N of 43 pix”^, significantly higher 
than our 10619 M31 spectra. The distribution of tem¬ 
plate spectral types is shown in Eigure[2 When referring 
to these spectra, we will use the term “template spectra” 
to distinguish them from the “science spectra.” 


2.4. AGB Identification 

We use position on the PHAT color-magnitude dia¬ 
grams (GMDs) identify stars at different evolutionary 
stages. The regions discussed in this section are illus¬ 
trated in Eigurej^ 

Stars on the red giant branch (RGB) were identi¬ 
fied using optical (F475IF and F814IF) Padova PAR- 
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F475W-F814W F110W-F160W 



a(A) 


Fig. 2. — Panel (a) shows an optical (F814VF vs. F475VF — F814VF) CMD of stars in the SPLASH sample, color-coded by evolutionary 
stage as determined from the photometry. Main sequence stars are plotted in blue, super giants in cyan, RGB stars in red, and AGB stars 
in pink. We show those stars identified as spectroscopically as carbon stars as yellow stars. Panel (b) shows the same stars on a NIR 
(F160VF vs. FllOVF — F160) GMD. Panel (c) shows a representative DEIMOS spectrum of each of the color-coded populations in the 
top panels. Spectra are normalized, smoothed by a Gaussian with cr ~ 3 pixels, and plotted with a vertical offset. We also mask the the 
telluric A band at 7600A. 


SEC1.2s isochrones (Bressan et al. 2012). The RGB 
boundaries, and most importantly tne tip of the red gi¬ 
ant branch (TRGB), were defined using isochrones with 
—2.1 < [M/H] < 0.6 and a fiducial age of 10 Gyr. We 
assumed a foreground reddening of Ay = 0.17 mag and a 
distance to M31 of 776 ± 18 kpc (jDalcanton et al.|2012|). 
For an example of the RGB footprint, see Figure 1 of 
G15. Super giants (SG) were identified as a distinct 
branch in F336VF - F475W vs. FllOVF - F160VF color- 
color space. All objects not classified as SG or RGB were 
classified based on optical color and magnitude. Those 
stars with F41bW — F8141T <1.2 are assumed to be 


main sequence (MS) stars, and those stars brighter than 
the TRGB with F41bW — F8141T >1.2 are assumed to 
be AGB stars. These designations are purely photomet¬ 
ric, with no reliance on spectral features. Our spectro¬ 
scopic sample is dominated by RGB stars (5680), with 
2489 MS stars and 1867 AGB stars. The remaining 583 
objects are classified as either SG or PNe. 

These designations are such that the bluest AGB stars 
may be bright supergiants and/or core helium-burning 
stars. We cannot distinguish between these objects and 
early-type AGB stars photometrically or spectroscopi- 
cally at our resolution (Melbourne et al. 2012 Dalcan- 
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Fig. 3. — Evaluation of the iF-statistic based carbon star identi¬ 
fication method. The solid histograms show the distribution of K 
values for our degraded templates with S/N> 2. For ease of visual¬ 
ization, we truncate the x-axis, and note that templates degraded 
to a relatively high S/N often have K values beyond these bounds. 


ton et al.]|2Q12 ). Alternatively, some may be the reddest 


(and relatively few) Milky Way stars tha t are observed 
at very red colors (| Williams et ah 2014). In addition, 
the distinction between KGB and AGB stars means our 
AGB sample is dominated by TP-AGB stars. 



2.5. Carbon Star Identification 

Our primary method of identifying carbon stars com¬ 
pares each science spectrum to our suite of 31 Milky Way 
template stars. To compute whether a science spectrum 
is best fit by a carbon template or a non-carbon tem¬ 
plate, we use a simple classification statistic K defined 
as: 


K = min (x^(6»jvc)) - min (x^(6»c)) (1) 

where Onc denotes the suite of non-carbon template 
stars, and Oc is the suite of carbon template stars. K 
is the difference between the statistic for the best fit¬ 
ting non-carbon template and the statistic for the best 
fitting carbon template. By this definition, those stars 
with K > are likely carbon stars. This metric is loosely 
based on the ratio of evidence, or Bayes Factor, with a 
minimization rather than a sum. 

We verify this method by degrading the S/N of each 
template to mimic the range of S/N present in our data, 
and then computing the degraded template’s K value by 
comparing it to the remaining 30 templates. This pro¬ 
cess is repeated until there are ^ 3000 degraded carbon 
templates and ^ 3000 degraded non-carbon templates. 
Figure shows the distribution of K for our degradation 
tests. They indicate that this method is ^ 88% accu¬ 
rate across the full range of S/N, and ^ 95% accurate 
for those stars with S/N > 2 pix“^. 

We apply our classification statistic to the 7903 stellar 
spectra with S/N pix“^ > 2, and visually inspect the 
spectra that return K > t). Specifically, we look for the 
GN features at ^ 7900A and C 2 features between 6100 
and 6600A, shown in Figure for the five carbon tem¬ 
plates and five randomly chosen carbon science spectra. 
We find a total of 94 carbon stars. There are 316 objects 


Fig. 4. — The left panel shows the five MW carbon stars used 
as templates for identification via iF-statistic. The right panel 
shows five M31 carbon stars chosen to demonstrate a range of S/N. 
Spectra are smoothed by a Gaussian with cr ~ 3 pixels, and the 
the telluric A band at 7600A is masked. The science spectra are 
labeled with their S/N pix“^. On both panels we have labeled the 
region where C 2 is prominent in blue (~ 6100 — 6600 A), and the 
region where CN is prominent in red ~ 7900 — 8200 A). These are 
the features we look for during the visual inspection that follows 
any automated carbon star detection method. 


whose spectra return TC > 0 but upon inspection do not 
have the carbon features highlighted in Figure ^ This is 
a higher false-positive rate than our tests predicted, be¬ 
cause our tests did not simulate the bad sky subtraction 
or untrustworthy radial velocity measurements common 
for a multi-slit spectrograph. 

To catch carbon stars that may have escaped identi¬ 
fication by the TC-statistic, we run a separate identifica¬ 
tion using synthetic narrow-band photometry modeled 
on the four band photometry system (FBPS). FBPS has 


Nowotny et al.||2003 

IBattinelli & Demers||2004a|ct |Wing| 

2007 

i Battinelli & ! 

Demers 2009b|). FBPS uses broad 


defined by the narrow band CN and TiO filters (centered 
at 8120.SA and 7778.4A, respectively) to separate C- and 
0-rich AGB stars. It is one of two preferred photomet¬ 
ric techniques for this separation (the other being NIR J 
and Ks photometry, whose reliability ha s recently been 
called into question; Menzies et al.|2015). 

To compute synthetic GIN and 'i'iO rnagnitudes from 
our spectra, we perform a first-order fiux calibration. We 
create a transmission curve by dividing our MW template 
spectrum of HD52005 with the fully flux-calibrated spec- 
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Fig. 5.— FBPS Color-color diagram (F475VF - F814VF vs. syn¬ 
thetic CN—TiO) used to identify potential carbon stars. The grey 
shaded region is defined by the bounds F475VF — F814VF > 2.0 
and CN-TiO> 0.3, and represents where carbon stars are likely 
to lie. Carbon stars identified by the K classification statistic are 
plotted as yellow stars, and carbon stars identified using this FBPS 
photometry (plus visual inspection of the spectra) are plotted as 
red stars. All other stars are plotted as black points. 


tmm of HD52005 from the X-shooter Spectral Library 
l2M4|, 
inape c 


(Chen et ah |2014|), and use this transmission curve to 
correct the shape of our spectra. We then weight the 
spectra by the CFHT/CFH12k CN and TiO through¬ 
put curves to generate synthetic photometryj^ We use 
F475bF — FSlAW as our broad band color. 

Because this method is computationally inexpensive, 
we run it on all stellar spectra to catch any carbon stars 
whose S/N excluded them from the LC-statistic method. 
The resulting FBPS diagram is shown in Figure We 
define the region of likely carbon stars using the position 
of carbon stars identified by the K classification statis¬ 
tic. The bounds are F475kF — FS14W > 2.0 and CN- 
TiO > 0.3, which is fully consistent with the bounds 
used throughout the literature. We visually inspect the 
spectra of each star in the bounded region, and find an 
additional 9 carbon stars. 

In total, we find 103 carbon stars in our sample of stars 
in the disk of M31. By using this combination of two 
automated-detection techniques plus visual inspection, 
we are confident that these 103 stars represent the vast 
majority of carbon stars in our spectroscopic sample. 

While the majority of these carbon stars lie predictably 
above the TRGB, four are found in other areas of the 
CMD (see Figure [^. Three have colors and magni¬ 
tudes matching the KGB population, and one lies in the 
MS. To verify that these stars are members of the M31 
disk rather than foreground MW carbon dwarfs, we com¬ 
pute their line-of-sight (LOS) velocities. We convert each 


® http://svo2.cab.inta-csic.es/svo/theory/fps3/ 


star’s redshift measurement to a velocity, then apply a 
telluric A-band correction and a correction for heliocen¬ 
tric motion on the date the spectrum w as obtained (for 
further details see [Dorman et al.|[^12[ and references 
therein). These four carbon stars have LOS velocities 
between —85 km/s and —300 km/s, making them likely 
M31 members. Their unusual colors and magnitudes 
are likely due to a mismatch between the SPLASH and 
PHAT catalogs. 

3. DEFINING M 


Our spectroscopic identification gives us an uncontami¬ 
nated sample of unambiguous carbon stars. There exists 
no equivalent spectroscopic selection that can separate 
M-type TP-AGB stars from the RGB. And while the 
G/M ratio has been computed across the Local Group, 
there is no standard photometric definition for M-stars. 
Even within M31, multiple criteria have been devised 
to isolate M-stars based on availabl e data. We examine 
three criteria b ased on the work by Boyer et al. (2013) 
(hereafter B13), Brewer et a l. ( 1995|) (hereafter B95), and 
Battinelli & Demers| (|2UU5|) (hereaiter BD05). Our goal 
is not to define a sin^e, “correct” method, but rather to 
investigate the differences in results returned by each of 
these methods. 


3.1. Boyer et al. 2013 

The first definition of M that we apply to our data is 
one modified from B13. As they use WFG3/IR medium- 
band filters to discriminate between G- and M-type TP- 
AGB stars, our dataset is certainly not identical. IR fil¬ 
ters are far more sensitive to dust-enshrouded stars than 
optical filters, and so the B13 sample will be substantially 
more complete. However, they make two major assump¬ 
tions that we can emulate: TP-AGB stars are brighter 
than TRGB, and all AGB stars not classified as G-stars 
are M-stars. Their definition of the TRGB also comes 
from PHAT photometry, so this identification of AGB 
stars matches our own. If we too take all non-carbon 
AGB stars to be M-type, we get a sample of 1605 stars 
from our original sample of 10619 spectra. 

Whenever this particular sample is used, it will be de¬ 
noted by a red diamond and will be referred to as the 
B13 method. 


3.2. Brewer et al. 1995 

B95 is one of the few large-scale AGB surveys of M31. 
The authors used FBPS in five regions located along the 
southern major axis. They define M by the following cri¬ 
teria: V - I > 1.8, GN-TiO < -0.2, / > 18.5, and 
^boi < —3.5 (where M^^/ = / + 0.3 + 0.34(E — I) — 
0.14(R —J)^ — (m —M)o). They assume a distance mod - 
ulus of 24.41, adopted from [Freedman & Madorej (|1990|). 

To replicate these criteria for our sample, we convert 
our FAlbW and F8141F magnitudes to V and I. To 
convert from FS14W to I we use the 97 M-g iant spec¬ 
tral templates presented by [Fluks et al.[ ([1994[) . We con¬ 
volve the Fluks et al. spectra with the HS1'/AGS F814VF 
bandpass and the GFIIT/GFH12k I bandpass. The off¬ 
set between F8141F and / is a function of color/spectral 
type. We fit a fourth-order polynomial to this relation¬ 
ship to determine a conversion function, getting residuals 
less than 0.002 magnitudes. To derive V — I we use the 
published transformations from Sirianni et al. (2005). 
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is used it will be denoted by a blue square and will be 
referred to as the B95 method. 
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F475W-F814W 


Fig. 6. — The results of applying the FBPS-based M-giant se¬ 
lection criteria of Brewer et. al 1995 to our spectroscopic sam¬ 
ple. The top panel includes all selection criteria; V — I > 1.8, 
CN-TiO < -0.2, I > 18.5, and Mboi < -3.5. Mboi = 
/ + 0.3 + 0.34(y - /) - 0.14(1/ -ly -{m-M)Q. The bottom panel 
includes all broad-band color and magnitude cuts but replaces the 
CN—TiO with the criterium that a star not be a spectroscopically- 
identified carbon star. AGB, RGB, and MS stars (denoted in grays) 
come from PHAT GMD-based identification. We adopt the cuts in 
the top panel for the B95 selection used in the rest of this paper. 


The application of B95 selection criteria is shown in 
Figure When we apply all selection criteria from B95, 
our final sample contains 843 M-stars (shown in blue in 
the top panel of Fig [^. The sample is relatively uncon¬ 
taminated, with only 2 stars having been identified by 
PHAT as RGB stars rather than AGB stars. However, it 
misses the bluest AGB stars (F4751F —F8141F < 3), the 
brightest AGB stars, and those stars close to the TRGB. 
The latter effect grows more pronounced at redder colors, 
and thus the more metal-rich populations. If we were to 
remove the GN—TiO criterium and replace it with the 
requirement that no st ar b e a spectroscopically identi¬ 
fied carbon star (see § |2.5[ ), the sample would increase 
to 1218 (shown in red m the bottom panel of Fig 
decreasing the resulting G/M ratio. This expands the 
sample to include many of the bluer AGB stars left out 
of the B95 sample. 

In the remainder of the paper, we will use the unmod¬ 
ified B95 selection of 843 stars. Whenever this sample 


3.3. Battinelli & Demers 2005 


A common way to express the G/M ratio in the Local 
Group is G/M0+; i.e. the M-type AGB stars counted 
are those with a spectral type of MO or later. This 
spectral-type selection is often defined photometrically 
using a color cut (e.g. R — I > 0.9). We can apply a 
color cut and determine spectral types from the spectra 
themselves. However, we must first note that the PHAT 
photometric observations and SPLASH spectroscopic ob¬ 
servations were taken several months (or in some cases 
years) apart. Since a significant fraction of TP-AGB 
stars are long-period variables of sizable amplitude, the 
measured PHAT colors may not refiect the colors of the 
star at the time the spectrum was taken. 

BD05 define M0+ to be those stars with R — / > 0.9 
and GN-TiO < 0. This broad-band color criterium 
matches that used by B95, as R — I = 0.9 corresponds 
to R — / = 1.8. To distinguish M-type AGB stars, they 
use a bolometric magnitude limit of Mboi < —3.5, where 
they adopt a slightly different definition of Mboi than 
B95: Mboi = / + 1.7103 - 2.2968(R - /) + 1.66464(R - 
J)^ — 0.43399(77 — — {m — M)o. BD05 adopt a dis¬ 

tance modulus (m — M)o = 24.41. To apply these same 
criteria we first transform F4751F and F814VF magni¬ 
tudes into / (using the method outlined in the previous 
section), and then de rive R using the transformations by 
(Sirianni et al.||2005|). The resulHng sample is shown in 
blue in the top panel of Figure It contains only 736 
stars, as the sample stops half a magnitude above the 
TRGB. This sample also omits the bluest AGB stars, 
but unlike B95, it does contain the brightest. As with 
the B95 criteria, if we were to replace the GN—TiO limit 
with the requirement that no star be a spectroscopically 
identified carbon star the sample would increase to 960 
stars (shown in red in the bottom panel of Fig [3, de- 
creasing the resulting G/M ratio. As before, this would 
extend the sample in the blue but would have no effect 
on the proximity of the sample to the TRGB. 

Since we have spectra in addition to photometry, we 
can also spectroscopically determine which stars have a 
subtype of MO or later. We use the TiO-based spectral 


parameters, and S 2 / 3 ,Sp defined in Fluks 


et al. (1994). These three parameters are calculated by 


integrating over the various TiO absorption features, and 
are monotonic functions of spectral type. Stars with sub- 
type MO or later have 5'i/2,5p < 0-517, < 0.754, 


and 5'2/3,Sp < 1-458. There is no spectroscopic distinc¬ 
tion between AGB and RGB stars at this resolution, 
so we again apply the bolometric magnitude limit from 
BD05. This returns a sample of 773 M-stars, shown in 
yellow on Figure The 736 M-stars identified by the 
photometric M0+ criteria are plotted in red. 

The sample of M-stars identified as M0+ by spectro¬ 
scopic criteria are almost identical to those identified by 
the photometric criteria. The top and right panels of 
Figure shows the marginal distributions of both the 
spectroscopic and photometric sample in color and mag¬ 
nitude. The color criteria are slightly less sensitive to 
fainter and bluer stars (F8141F > 18). 

In the remainder of the paper, we will use the unmod- 
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Fig. 7.— The results of applying the FBPS-based M-giant selec¬ 
tion criteria of Battinelli & Demers 2005 to our spectroscopic sam¬ 
ple. The top panel includes all BD05 selection criteria; R — I > 0.9, 
CN-TiO < 0, and Mboi < -3.5. Here, Mboi = I + 1.7103 - 
2.2968(ii - /) + 1.66464(J? - 7)2 - 0.43399(7J - if - (m - M)o. 
The bottom panel replaces the CN—TiO criteria with the require¬ 
ment that a star not be a spectroscopically identified carbon star. 
AGB, RGB, and MS stars (denoted in grays) come from PHAT 
GMD-based identification. We adopt the cuts in the top panel for 
the B05 selection used in the rest of this paper. 


ified BD05 photometric selection of 736 stars. Whenever 
this sample is used it will be denoted by a green circle 
and will be referred to as the BD05 method. 


4. G/M RATIO AGROSS M31 

The broad range of environments in M31 mean that we 
can compute the C/M ratio as a function of a variety of 
properties. In this section we look at th e C/ M ratio as a 
fu nctio n of galactoc entr ic radius (i?g, § 4.1), metallicity 
(§ 4.2), and SFH (§ |4.3[) in spatial bins. 


Ae create the spatial bins used throughout this section 
beginning with equal-area bins on a flat circle with a 
radius of 20 kpc (roughly the outer limit of our data). 
We then incline the bins to match M31’s inclination of 


74° (jBarmby et al.j2006|) and tilt the bins to match M31’s 
position angle of 5m (determined empirically). The flnal 
positions of these bins are shown along with the positions 


of our C- and M-stars in Figure 

4.1. ^45 a function of galactocentric radius 

Using the bins shown in Figure we determine the 
change in the C/M ratio as a function of radius. We 
compute the center radius of each bin by taking the av¬ 
erage Rg of all PHAT stars within that bin. When com¬ 
puting this average, we leave out PHAT “Brick 1.” Brick 
1 covers the extremely crowded bulge region, extending 
only ^2.7 kpc from the center of M31, where we have 
no spectroscopic coverage. Leaving those stars out of our 
average ensures that the resulting distances represent the 
regions in which our AGB stars lie. The resulting C/M 
ratios for each of the selection criteria discussed in ^are 
shown in Figure pU) 

There is a clear positive trend. A weighted least 
squares fit to the eight bins along major axis bins de¬ 
termines the following gradients: 

log(C/MBi 3 ) = (0.06±0.016)xi?g-(1.92±0.22) 
log(C/MB 95 ) = (0.05±0.016)xi?g-(1.56±0.23) 
log(C/MBD 05 ) = (0.06±0.017)xi?g-(1.53±0.23) 

This trend is likely to be driven not by the radius itself, 
but by other quantities that vary with radius. We con¬ 
sider correlations with these quantities in the following 
two sections. 


4.2. As a function of metallicity 

We are in a unique position to compute the C/M ra¬ 
tio as a function of metallicity calculated using different 
techniques and for different tracer populations. In this 
section we focus on gas-phase metallicity, as determined 
from HI I regions, and on photometrically derived stellar 
metallicity. 



X 
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Fig. 8. — The result of a spectrum-based M-giant selection cri¬ 
terion applied to our sample. As in Figures and AGB, RGB, 
and MS stars are in greys. Yellow crosses represent those stars 
with spectral types of M O or later, as determined by the indices 
from |Fluks et al.j ( |1994| > 
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Fig. 9. — Positions of all AGB stars in the SPLASH sample (pink 
points), with carbon stars (yellow stars), overl aid on a GAL EX 
image of M31. The spatial bins used throughout §4.1| §4. 2 1 and §4.3| 
are marked with black lines. Apparent structure m the positions 
of the AGB stars is due to the SPLASH selection function, and is 
not a true physical effect. 


4.2.1. Gas-phase metallicity 

We first look at the gas-phase metallicities from 
Sanders et al. (2012). The authors measure the oxy- 
gen abundance ( |0/H]) in HIT rejgious u sing strong-line 
diagnostics from Zaritsky et al . (199 4|), Kewley fc Do- 
ita (|2QQ2|), jlNagao et al. (|2UU6p, and |r ilyugin ^ T'huan 
2005 r They hnd the hrst three m ethods very consis¬ 
tent, p ut note that the method from |P ilyugin fc Thuan 
(2005) has a median offset ^ 5x that of the others. In 
l ight o f this, when compiling the data from [Sanders et aT 
(2012|) we disregard the measurements from Pilyugin & 


.'huan (2005) and take a weighted average of the rest, 
.'he positions of the HI I regions are shown with respect 
to our data in Figure |TTj 

We compute the C/M ratio as a function of gas-phase 
[0/H] using the bins shown in Figure 1^ The [0/H] 
value for each bin is defined as the [0/H] value of the 
HI I region closest to the average position of all PH AT 
stars within the bin. We chose this method rather than 


an interpolation because Sanders et al. (2012) find that 
the [O/H] distribution is quite clumpy. While there is a 
global metallicity gradient of —0.0195T0.0055 dex kpc“^, 


^ B13 method 

M B95 method 
0 BD05 method 





■^■^4 6 8 10 12 14 16 18 

R (kpc) 


Fig. 10.— Log(G/M) ratio in bins of Rg. The radii are those 
of the mean stellar distance in each bin. Red diamonds represent 
points calculated using the B13-based definition of M, blue squares 
represent the B95-based definition of M, and green diamonds rep¬ 
resent the BD05-based definition of M. For ease of visualization, 
B95 and BD05 points are artificially offset from the B13 points 
by 0.2 and -0.2 kpc, respectively. Bins along the major axis are 
shown as solid points and bins along the minor axis are shown as 
transparent points. 


HI I regions even ^ 2 kpc apart may have dramatically 
different abundances. 

Panel (a) of Figure 12 show the resulting relationship 
between log(C/M) and [0/H]. Panel (b) shows the same 
relationship for only those bins along the major axis. 
There is a negative gradient, with the C/M ratio decreas¬ 
ing with increasing metallicity. We find the weighted 
best-fit to the major axis bins to be: 


log(C/MBi 3 ) = (-0.96±0.28)x[0/H]+ (7.31T2.48) 
log(C/MB 95 ) = (-0.95±0.29)x[O/H] + (7.53±2.55) 
log(C/MBD 05 ) = (-0.90±0.29)x[0/H] + (7.10±2.58) 


HI I regions probe the metallicity of present-day star 
forming regions and young stars, not the relevant frac- 



Fig. 11.— Positions of HII regions (circles) and PNe (triangles) 
from Sanders et al. 2012, with respect to our spatial bins (black 
lines). Points are colored by their oxygen abundance. 
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[0/H] [0/H] 




[M/H] [M/H] 


Fig. 12.— Log(C/M) vs various properties of the M31 disk. In each panel, red diamonds are calculated using the B13-based definition of 
M, blue squares represent the B95-based definition of M, and green diamonds represent the BD05-based definition of M. Panel (a) shows 
the C/M ratio vs. approximate [O/H] value at the center of each spatial bin. Panel (b) shows the C/M ratio vs approximate [O/H] value 
at the center of bins along the major axis only. Panel (c) shows the C/M ratio vs. mean photometric metallicity for all bins. Panel (d) 
shows the C/M ratio vs. mean photometric metallicity for bins along the major axis only. 


tion of AGB stars that formed several Gyr ago. As a 
result, correlations between [0/H] in HII regions and 
the G/M ratio tell us more about the metallicity evolu¬ 
tion of the galaxy than the formation and evolution of 
the AGB stars themselves. An alternative source of gas- 
phase oxygen abundances are planetary nebulae (PNe), 
w hose ages are much c loser to those of AGB stars. 


Sanders et al. (|2012) compute the oxygen abundance 
of P^e in the disk of M31. Unfortunately, there are not 
enough PNe in the vicinity of our data for us to look at 
the C/M ratio as a function of thei r [O/H] measur e ments 
(see Figure pT]). Recen t work by Kwitter et al. (2012) 
and |Balick etal. (2013) compute oxygen abundances lor 


PNe in the outer disk of M31 {Rg > 18 kpc). Upcoming 
spectrosc opic follow-up of th e many PNe in the PH AT 
footprint (Veyette et al. |2014| ) will provide reliable [O/H] 
measurements against which we can compare the C/M 
ratio. 


4.2.2. Stellar metallicity 

In addition to a gas-phase metallicity, we can study the 
C/M ratio as a function of stellar metallicity. We use 
the photometric metallicity estimates derived for RGB 
stars in the PH AT fields by G15. While the photomet¬ 
ric metallicity of the RGB population may not exactly 
match the metallicity of the AGB population, it is a far 
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^AGB^RGB ^AGB^RGB 




Mean SFR over past 400 Myr (M yr Mean SFR over past 400 Myr (M , yr 


Fig. 12 (CONT.). — Panel (e) shows the C/M ratio vs. our age proxy, ^ agbRGB ^ for all spatial bins. Panel (f) shows the C/M ratio 
against bins along the major axis only. Panel (g) shows the C/M ratio vs. the mean SFR over the past 400 Myr for all spatial bins. Panel 
(h) shows the C/M ratio against the mean recent SFR for bins along the major axis only. 


better tracer than HI I regions. In addition, we find that 
the combination of many more stars and more robust 
models makes photometric metallicity estimates derived 
from the RGB stars much more reliable than those we 
can derive using the spectra of our AGB stars. 

G15 determine stellar metallicities by interpolating 
bias and completeness corr ected photometry ont o Padova 
PARSEG1.2s isochrones (Bressan et al. 2012), with a 
fixed age of 4 Gyr and a metallicity range —2.18 < 
[M/H] < 0.6 (0.0001 < Z < 0.06). To avoid being heav¬ 
ily biased by very crowded regions, they restrict their 
interpolation to those stars with F814IT < 23. In the 
outer regions of the disk this cut corresponds to 100% 
completeness, while in innermost regions it marks ^ 50% 
completeness. Within the spatial bins shown in Figure 


we compute the mean metallicity of all RGB stars satis¬ 
fying the aforementioned magnitude and metallicity cri¬ 
teria. We leave out the metallicities of stars in PHAT 
Brick 1, where we have no spectra, so as not to artifi¬ 
cially boost the average metallicity of the innermost bins 
beyond what is representative of our spectral data. 

Figure panel (c), shows log(G/M) as a function of 
mean metallicity for all spatial bins. Panel (d) shows the 
relationship for only those bins along the major axis. We 
see a statistically significant negative gradient in G/M 
with photometric metallicity. A weighted best-fit to the 
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spatial bins along the major axis returns 
log(C/MBi 3 ) = (-2.84±0.75)x[M/H]-(1.21±0.06) 
log(C/MB 95 ) = (-2.63±0.78)x[M/H]-(0.95±0.07) 
log(C/MBD 05 ) = (-2.71±0.79)x[M/H]-(0.89±0.07) 

G15 note that the radial metallicity gradient computed 
from photometric RGB metallicities depends on the un¬ 
derlying age distribution, and provide median metallicity 
measurements for several fiducial ages. To test the de¬ 
pendence of the G/M ratio on age assumptions, we com¬ 
pute the gradient in log(G/M) with photometric metal¬ 
licity for fiducial ages of 6 Gyr and 8 Gyr. For the B13 
method, we find slopes of—2.78 ±0.73 and —2.73 ±0.72, 
respectively, which is perfectly consistent with the slope 
determined for a fiducial age of 4 Gyr. 

In addition to depending on the underlying age dis¬ 
tribution, the metallicity determinations from G15 are 
affected by the choice of isochrones, low level dust ex¬ 
tinction throughout the survey region, and any uncor¬ 
rected photometric bias and completeness. As a result, 
the absolute metallicity measurements have large uncer¬ 
tainties. 


4.3. As a function of SFH 

In addition to being a function of metallicity, the G/M 
ratio should depend on a galaxy’s SFH. For a single-burst 
stellar population, the G/M ratio begins to increase from 
zero once turn-off masses reach 3-4 Mq. For ages older 
than a few Gyr, carbon stars are no longer produced and 
the G/M ratio falls again to zero. The exact age at which 
the G/M ratio peaks is determined by star formation rate 


(SFR ) and metallicity (Gioni et al. 2006 Marigo et al. 

In a galaxy, this behavior is convolved with the 
S4'H and smoothed. Using PHAT data, we can evaluate 
the G/M ratio as a function of population age and SFR. 

While the exact ages of the AGB stars in M31 are 
difficult to compute, we can use the PHAT data to define 
a rough proxy. Again, we use the spatial bins shown 
in Figure We define our proxy to be the number of 
TP-AGB stars in a given bin divided by the number of 
RGB stars The value of this quantity as 

an age proxy relies on the fact that the mean progenitor 
mass of TP-AGB stars is larger, and hence younger, than 
that of RGB stars. A higher number ratio of TP-AGB 
stars to RGB stars implies a younger average population. 

In order to be sensitive to the relative numbers of AGB 
and RGB stars independent of photometric complete¬ 
ness, we set a lower magnitude limit of F8I4IU = 23. 
We also leave out the data from PHAT Brick I, so as 
not to be biased by the crowded bulge in which we 
have no spectra. With the crowded regions inward of 
Rg ^ 4 kpc left out, our magnitude limit represents 
100% completeness in our bins. We also set a color limit 
of FAlbW — F8I4IU > 2 to mitiga te co ntamination by 
younger helium burning stars (see § 2.4). 

A plot of our age proxy as a function of Rg is shown in 
Figure As this verification does not require SPLASH 
data, whose smaller numbers limit how small we can bin. 
Figure is made with polar bins created by splitting 
M31 into 31 equal-area bins rather than the thirteen 
shown in Figure As expected in the regime of hier¬ 
archical galaxy formation, we see that populations are 
younger as you move farther out in the disk. 



Fig. 13. — Age proxy agb/^rgb) as a function of Rg in kpc. 
Error bars represent Icr Poisson uncertainties. 


Figure 12, panel (e), shows log(G/M) as a function 
of this ^e proxy computed in all spatial bins shown in 
Figure]^. Panel (f) shows the same for bins along the 
major axis. The weighted best-fit lines to bins along the 
major axis are: 


log(G/MBi3) = ( 22 . 5 ± 6 . 5 )X 
log(G/MB95) = (I8.8±6.6)x 
log(G/MBD05) = (19.I±6.6)x 


Nagb 

Nrgb 

Nagb 

Nrgb 

Nagb 

Nrgb 


(2.46±0.40) 

(I.98±0.40) 

(I.93±0.40) 


We next look at th e G/M ratio as a fu nction of recent 
SFR, using data from|Lewis et al. (20I5|) (hereaf ter LI5). 
LI5 a pply the GMD fitting program MATGH ( [Dolphin 


2002) to PHAT data in over 9000 regions in the disk 


of M31, and compute SFRs and cumulative stellar mass 
formation. To minimize the effect of model uncertain¬ 
ties their analysis focuses solely on the main sequence, 
and so probes the SFH of M3I over the last 400 Myr. 
While this timescale does not match the period during 
which the bulk of the AGB stars formed, both models 
and observations predict TP-AGB stars in populations 
as young as 100 Myr (e.g. [Frogel et al.||I990|). 

We take each of the 9000 regions investigated by LI5 
and group them by the spatial bins shown in Figure to 
compute a mean SFR. Log(G/M) as a function of mean 
recent SFR is shown in Figure for all bins (panel g) 
and for bins along the major axis (panel h). There is 
considerably more scatter in these relationships than we 
saw for gas-phase oxygen abundance, stellar metallicity, 
or age proxy. However a weighted best-fit line to ma¬ 
jor axis bins indicates a statistically significant negative 
correlation. 


log(G/MBi3) = (-238±95)xSFR-(0.9I±0.09) 
log(G/MB95) = (-229±98)xSFR-(0.67±0.09) 
log(G/MBD05) = (-24I±99)xSFR- (0.59±0.09) 


5. DISCUSSION 

5.1. The C/M ratio vs. environment 
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TABLE 1 

Coefficients of Multiple Regression Models 



[O/H] 

[M/H] 


SFR 

R 

Model 1 


-2.25 ±0.95 

9.61 ±7.7 



Model 2 

-0.59 ±0.24 

-0.99 ±0.86 

11.4 ±4.9 



Model 3 


-1.79 ± 1.48 

12.5 ±9.8 

-67.6 ± 149 


Model 4 


-1.72 ± 1.54 



0.04 ±0.03 

Model 5 

-0.64 ±0.17 




0.05 ±0.01 

Model 5 


-2.38 ± 1.92 

10.6 ±16.2 


-0.01 ±0.07 


Note. — Rows separate individual regression models, and columns represent which 
parameters were used in each model (gas-phase oxygen abundance [O/H], stellar metal- 
licity [M/H], age ^ agbRGB : mean SFR over the past 400 Myr, or galactocentric 
radius Rg). A blank entry indicates that the parameter was not used in that particular 
regression model. 


We find statistically significant trends between 
log(C/M) and each of the five properties we investi¬ 
gate: i?g, present-day gas-phase [O/H], stellar [M/H], 
age proxy, and recent SFR. In each case, the gradient is 
made clearer by restricting the analysis to the eight spa¬ 
tial bins along the major axis. This is likely due to the 
increased uncertainty in the deprojection along the mi¬ 
nor axis rather than an azimuthal change in the behavior 
of the C/M ratio. 

The clearest trend we see is between log(C/M) and ra¬ 
dius. Our slopes 0.06) are perfectly in line with the 
gradients observed by B95 (excluding thei r measurement 
at ^ 30 kpc) and Battinelli et al. (2003) along the op¬ 
posite axis of the M31 disk. While there is no reason 
the C/M ratio should depend on geometry, this correla¬ 
tion points to a significant relationship between the C/M 
ratio and properties that themselves vary with radius. 

B13 measure the C/M ratio in the inner disk of M31 
{Rg = 2 kpc), and find a C/M ratio much lower than 
predicted given the metallicity at that radius. They in¬ 
terpret this as an indication of a metallicity ceiling above 
which carbon stars do not form. If we extrapolate our 
radial gradient to Rg = 2 kpc, we calculate log(C/M) = 
— 1.81 ± 0.22. This is considerably higher than the 
value measured by BI3 (log(C/M)= —3.48lo'.oi)* 
paucity of carbon stars observed in the inner disk by 
BI3 was due to a gradual fall off in carbon star produc¬ 
tion with increasing metallicity, then our radial gradient 
should extrapolate to their measurement. The fact that 
the C/M ratio measured by BI3 is so low is consistent 
with the suggestion that there is a hard metallicity limit 
above which carbon stars do not form. 

While stellar oxygen abundance is the major driver of 
whether a TP-AGB star can become a carbon star, it 
is interesting that the [O/H] measurements from HII re¬ 
gions in the disk of M31 show the expected trend when 
they themselves do not probe the metallicity of the AGB 
stars. One possible explanation for the observed gradi¬ 
ent is that the metallicities of present-day star forming 
regions are correlated with the metallicities of star form¬ 
ing regions a few Gyr ago. This would indicate that the 
metallicity gradient in M3I is long lived. Alternatively, 
the G/M ratio and present-day gas-phase [O/H] may be 
independent, and the observed gradient is an artifact of 
their shared relationship with radius. 

The clear gradients we observe in log(G/M) with stel¬ 
lar metallicity and age are predicted by models and pre- 


served relationship between log(G/M) and age indicates 
that young populations produce more M- than G-type 
TP-AGB stars, which could be explained by hot bottom 
burning operating on the most massive TP-AGB stars. 
The fact that the G/M ratio increases monotonically with 
our age proxy indicates that the average age of the pop¬ 
ulations in all our bins are older than a few Gyr. If this 
were not the case, then we would see the G/M ratio turn 
over and begin to decrease at the bins dominated by very 
young populations, where stars evolving off the main se¬ 
quence do not become AGB stars or remain M-type due 
to hot bottom burning. Average ages consistently greater 
than a few Gyr is consistent with recent work showing a 
disk-wide b urst of star formation 2-4 Gyr ago ([Willian^ 


et al.|j20I5), and the age calculations by Dorman et al. 
(|2UI5| and GI5. 

Our derived trend with metallicity is fully consistent 
with the relationship established acros s the Local Group 
(log(G/M)- (-2.12 ±0.04) X [Fe/H]; |Gioni||2009D. The 
trend observed in Local Group satellites appears to ex¬ 
tend smoothly to the metal-rich regime of M31, irrespec¬ 
tive of their drastically different SFHs. This implies that 
the G/M ratio is strongly determined by the same popu¬ 
lation of intermediate-age stars in all environments, and 
depends only on stellar properties (metallicity and age). 

This interpretation of the log(G/M)-metalhcity rela¬ 
tionship would imply that our observed gradient in the 
G/M ratio with recent SFR is driven by correlations be¬ 
tween SFR, age and metallicity rather than the impact of 
SFR itself on AGB evolution. The alternative explana¬ 
tion for that trend is that star formation in M31 is long 
lived, and the SFR of the past 400 Myr correlates with 
the SFR a few Gyr ago. This correlation coul d occur if 
the global b urst of star formation 2-4 Gyr ago (Williams 
et al [2015 ) was spatially heterogeneous, and the regions 


vious observations (e.g. Mouhcine fc Lancon||2003 


et al.|201Q Held et al.|20lO| Marigo et al.|20l'3 ). Our 


Feast 


of the disk producing the most stars 2-4 Gyr ago also 
have the highest recent SFR. 

To investigate whether we can decouple the effects of 
metallicity and age on the G/M ratio and determine the 
dominant factor, we construct a series of multiple regres¬ 
sion models. Beginning with a simple linear relationship 
between log(G/M) and any property of the disk, analysis 
of variance finds that the model improvement from the 
addition of a second (or third) variable is insignificant. 
The coefficients of a sample of these models are shown 
in Table none are non-zero with high confidence. This 
is indicative of collinearity, and that age and metallicity 
are so tightly coupled that we cannot claim that statis¬ 
tically one is more important to the G/M ratio than the 
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other. 


5.2. Possible effects of dust 

Both ISM extinction within M31 and circumstellar 
dust may affect the measurement of the C/M ratio, and 
ISM extinction may affect our metallicity, age proxy and 
SFR measurements. 

The C/M ratio will be affected by dust if it causes car¬ 
bon stars or M-stars to be preferentially excluded from 
the SPLASH sample. Since we identify carbon stars 
spectroscopically, veiling of the prominent CN and C 2 
features would make the stars unidentifiable. This is not 
the case for M-stars, as we identify them photometrically. 
However, recent models show that even the dustiest car¬ 
bon stars do not have their r\j 79 OOA CN feature com¬ 
pletely veiled (Aringer et. al, private communication). It 
is far more likely that carbon stars escape our selection 
by being too dust-reddened to be visible in the optical, 
and are thus not targeted for spectroscopy. 

To investigate the effects of interste l lar du st, we use the 
M3I dust maps from Dalcanton et al. (2015) to determine 
the average extinction [Ay) at the position of each star. 
We compare the distributions of Ay at the locations of 
the carbon stars and M-stars via a Kolmogorov-Smirnov 
(KS) test, and find no significant difference. Interstellar 
dust is thus not likely to have a significant effect on our 
selection. 

Circumstellar dust is known to have a major i mpact 
on the comp leteness of optical AGB surveys (e.g. Boyer 


et al. 2011). Correcting for AGB stars reddened out 
of the optical by circumstellar dust would steepen our 
measured age and metallicity gradients. In the case of 
metallicity, this is because while dust production in car¬ 
bon stars stays constant as me tallicity increases, dust 
production in M-stars goes up (Sloan et al. 2008). As 
metallicity increases the number of M-stars reddened 
from our sample will increase faster than the number 
of carbon stars reddened from our sample. The gradi¬ 
ent with age proxy will steepen because carbon stars are 
typically more massive than M-stars, and produce more 
dust. Thus as the population gets younger, and has more 
massive stars climbing the AGB, the fraction of carbon 
stars invisible in the optical will be increasingly larger 
than the fraction of M-stars invisible in the optical. 

Because both age and metallicity decrease with radius 
in M31, both effects discussed above will contribute to 
our measuring a shallower gradient than reality. 

By analyzing the Ay and the fraction of reddened stars 
ifred) GI5 determine that their median metallicity mea¬ 
surements are generally not dependent on dust. They do 
find that the high metallicities found for stars along the 
north-west edge of the PHAT footprint track regions of 
higher dust extinction, however our spatial bins are large 
enough that these regions make up a small fraction of the 
stars do not substantially bias our median metallicities. 

To evaluate the effect of dust on our age proxy, we 
first compare the distribution of extinction values in the 
vicinity of the AGB and R GB stars of the PHAT sa mple 
using the dust map from jDalcanton et al.| (|20I5|). A 
KS-test indicates that we can reject the hypothesis that 
these populations are drawn from the same distribution 
of Ay. AGB stars are found in slightly dustier regions. 
However, we find no trend between age proxy and Ay 
(r 2 = 0.09), which indicates that this does not heavily 


impact our results. 

LI5 have taken dust into account when determining 
SFRs, so we will not further analyze the effect of dust on 
these measurements. 


5.3. Comparing definitions of M 

Throughout Section we see that the definition of M 
in the G/M ratio has a significant impact on the intercept 
of gradients but a minimal impact on the slope. Gonsis- 
tently, the B95 and BD05 methods return very similar 
results, despite having rather different effective TRGBs. 
The B13 method, on the other hand, always counts more 
M giants and returns a lower G/M ratio. This result in¬ 
dicates that care must be taken when comparing values 
within the literature, as inhomogeneous methods lead 
to different absolute values of log(G/M). The fact that 
slopes are consistent between different selection methods 
while intercepts are not tells us that the population of 
stars not being included - the bluer “AGB” stars - are 
distributed evenly across the disk. As we see gradients in 
log(G/M) with both age proxy and metallicity, this blue 
population is likely not the youngest or most metal poor 
AGB stars, but is instead either foreground contamina¬ 
tion or supergiants. 

The major result from B13, that the G/M ratio cal¬ 
culated ^ 2 kpc from the center of M31 is highly dis¬ 
crepant with values calculated farther out in the disk, 
stems from a comparison between the G/M ratios calcu¬ 
lated by the BI3 and B95 methods. B13 thoroughly in¬ 
vestigate whether this discrepancy stems from an under¬ 
estimation of the number of G-stars, but do not consider 
whether they are overestimating the number of M-stars 
(relative to the studies to which they compare). We can 
estimate how much the BI3 data point would change if 
the B95 criteria were applied. In our data, using the 
BI3 definition rather than the B95 definition amounts to 
a difference of ^ 800 M stars, or a difference of ^ 0.35 in 
log(G/M). While this moves the BI3 point closer to the 
trend fit to B95 data, it is still significantly discrepant (off 
by ^ I in log (G/M)). However, this rough approximation 
does not take into account the fact that B13 use medium- 
band NIR photometry to count M-stars while we use op¬ 
tical photometry. We cannot reproduce the effectiveness 
of medium-band NIR photometry, which is centered on 
distinctive absorption features present in either G or M 
stars, with PHAT wide-band NIR photometry. 


6. CONCLUSION 

We have computed the G/M ratio in the disk of 
M31 using spectroscopic and photometric data from 
the SPLASH and PHAT surveys. An uncontaminated 
sample of carbon stars was identified using moderate- 
resolution optical spectra. M-stars were identified pho- 
tometric ally in three differe n t ways, followi n g met hods 
used by Boyer et al . (1201 3), Brewer et al. (1995) and 
attinelli Demers 112005 ). We have calculated the G/M 
ratio as a function of galactocentric radius, present-day 
gas-phase metallicity, stellar metallicity, age (via proxy 
NAGB/Ni?GB), and mean SFR over the past 400 Myr. 

From this, we conclude: 


• The definition of “M” has a minimal effect on the 
slope of a relationship, but a substantial effect on 
the fiducial value. This adds to the body of evi- 
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dence that stresses the need for homogenous sam¬ 
ples when studying the C/M ratio as a function of 
other parameters. 


• There are statistically significant correlations be¬ 
tween log(C/M) and i?g, present-day gas-phase 
[0/H], stellar [M/H], age proxy, and recent SFR. 


We reproduce the relationship between log(C /M) 
and stellar metallicity stated in Cioni (2009|), de¬ 
spite working in a substantially more metal-rich en¬ 
vironment with a drastically different SFH. 


• Age and metallicity are too closely connected to 
state that one is more important to the C/M ratio 
in M31 than the other. 
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